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Abstract 

Charmonium states are considered important signatures of the strongly interacting medium cre- 
ated in heavy-ion collisions. In the ALICE experiment, these probes can be investigated in the 
pfpT decay channel, in the forward rapidity region (2.5< y <4) down to zero transverse mo- 
mentum. Results on charmonia production in Pb-Pb collisions at y/s NN -2.16 TeV are presented. 
The centrality and transverse momentum dependence of the inclusive J /iff nuclear modification 
factor are shown and compared with theoretical models. Finally, first ALICE results on the i/'(2S) 
production in Pb-Pb collisions are also discussed. 



The ALICE experiment at the Large Hadron Collider (LHC) is designed to study the forma- 
tion of the strongly interacting matter in high energy heavy-ion collisions. Among the probes 
of the expected phase transition between hadronic and deconfined matter, a relevant role is 
played by quarkonium states. In particular, according to the color-screening model [1], the in- 
medium dissociation probability of such states should provide an estimate of the initial temper- 
ature reached in the collisions. Studies performed in the last twenty years at the SPS and RHIC 
facilities, at ysAw=17 and 200 GeV respectively, have, indeed, shown a reduction of the J/i/f 
production yield beyond the expectations due to the cold nuclear matter effects (i.e shadowing 
and nuclear absorption). In spite of the very different center of mass energy, the amount of sup- 
pression observed by SPS and RHIC experiments is similar. Furthermore at RHIC a stronger J/iff 
suppression has been measured at forward with respect to mid-rapidity. These observations sug- 
gest the existence of an additional J/iff production mechanism, which sets in when higher y/s^N 
are reached. This mechanism, based on the combination of initially uncorrelated c and c pairs 
EES!, can counteract the quarkonium suppression in the QGP Therefore, the measurement of 
charmonium production is especially promising at the LHC, where the high-energy density of 
the medium and the large number of cc pairs produced in central Pb-Pb collisions should help to 
disentangle suppression and regeneration scenarios. 

The ALICE experiment [5 1 studies quarkonium production in the pfpT decay channel in the 
rapidity region 2.5 < y < 4 and in the <? + e~ decay channel at \y\ < 0.9 J6). In this paper, we present 
quarkonium measurement at forward y, while results obtained at mid-rapidity are discussed in 
0. 

Muons are identified and tracked in a Muon Spectrometer equipped with a dipole magnet, 
a set of absorbers, five tracking chambers and a trigger system. The pixel layers of the Inner 
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Tracking System allow the vertex determination, while forward detectors such as the VZERO 
scintillators are used for triggering purposes. The VZERO is also used to determine the centrality 
of the collisions, through a fit, based on a Glauber model, to its signal amplitude [8 |. Results 
presented in the following are based on a sample of dimuon triggered events corresponding to 
an integrated luminosity L ~ 70 /Ar 1 , collected during the 2011 Pb-Pb data taking. A clean 
data sample is obtained applying quality cuts to the muon tracks. These cuts eliminate mainly 
tracks hitting the edges of the spectrometer's acceptance or crossing the thicker part of the beam 
shield. Furthermore, to remove hadrons produced in the front absorber, tracks reconstructed in 
the tracking chambers are required to match hits in the trigger planes. 

The in-medium modification of the J/i// production is quantified through the nuclear mod- 
ification factor (Raa), defined as the ratio of the J/tff yield measured in Pb-Pb collisions and 
the expected yield obtained scaling the pp J/i// production cross section by the number of bi- 
nary nucleon-nucleon collisions. The ]/iff yield is extracted by fitting the opposite sign invariant 
mass spectrum with an extended Crystal-Ball function (CB2), which allows the inclusion of non- 
gaussian tails on both sides of the ]/tf/ pole. The mass position and the width of the CB2 are kept 
as free parameters in the fit, while the tails are tuned on a MC where a }/if/ signal is embedded into 
real events. The background contribution under the J/i/t is described by a gaussian function with a 
mass-dependent width. Alternatively, an event mixing procedure has been applied, to subtract the 
background contribution, before fitting the signal. The raw i/ift yield is determined as the average 
of the results obtained with the two approaches, including also some modifications of the fitting 
procedure (e.g. different CB2 sets of tails or alternative fitting functions). The corresponding 
systematic uncertainties on the signal extraction are defined as the r.m.s. of these results. Details 
on the analysis are given in [9|. The total number of }/if/ in the kinematic region 0< pj <8 GeV/c, 
2.5< y <4 and in the centrality range 0-90% amounts to ~40000. The high statistics collected in 
2011 allows, therefore, a differential study of the J/i/f production yield, as a function of central- 
ity, y and pj. The procedure described above is applied in each kinematic bin under study. To 
evaluate the Raa the J/i/f yield is divided by the acceptance x efficiency (A x e), computed by 
embedding generated J/if/ particles into real events. The average A x e is ~14%, with an 8% de- 
crease from peripheral to central collisions. Finally, the i/ij/ yield measured in Pb-Pb collisions, 
in each kinematic bin, is normalized to the corresponding inclusive i/ip cross section measured 
in pp collisions at the same energy iflOl . Systematic uncertainties on the Raa depend on the 
kinematic range under study. The main source is due to the pp reference (~9%), while other con- 
tributions, related, for example, to the choice of the MC inputs, to the uncertainty on the trigger, 
tracking and matching efficiency and to the signal extraction amount to less than 6-7% each. The 
inclusive Raa integrated over centrality, p T and y is A 90% = 0.497 ± 0.006(sfaf) + 0.078(sysf), 
exhibiting a clear i/ifr suppression. In our y and pj domain, the contribution from beauty hadron 
feed-down to the inclusive J/if/ yield amounts to ~10%, having a negligible effect on the Raa 
measurement. The centrality dependence of the Raa, integrated over pj and y is shown in Fig[T] 
(left). The pattern observed by ALICE presents a weaker centrality dependence and a smaller 
suppression for central collisions with respect to PHENIX results IfTTI . suggesting a different 
interplay, at the two energies, of suppression and regeneration mechanisms. Models including 
a large fraction of J/if/ produced from pair (re)combination lfT2l [T3l [l4ll or all J/i/r produced at 
hadronization [3| can reasonably describe the measured Raa- The (re)combination contribution 
is expected to be dominant, especially in central collisions, at low p T , while it becomes negligible 
as the J/if/ p T increases. This behaviour can be investigated by further studying the Raa Pi de- 
pendence in centrality bins, as shown in Fig [T] (right). While for the peripheral bin (40-90%) the 
pj dependence is negligible, in the most central bin (0-20%) the suppression increases by ~60% 
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Figure 1: Left: the J/i// Raa is shown as a function of Wpart and compared to the PHENIX result. Statistical errors are 
shown as lines, while systematic uncertainties are shown as boxes around the points. Right: the Raa Pt dependence 
is shown for two centrality classes and compared to a theoretical model 1 12']. In this case, systematic uncertainties are 
split between uncorrelated and partially correlated components, shown, respectively, as boxes and brackets around the 
points. The theoretical bands represent the uncertainty on the calculation due to different shadowing assumptions. Global 
systematic uncertainties are listed in the legend. 



moving towards low p T . Models including a /? T -dependent contribution from (re)combination 
which amounts, at low pj, to 60% in central and 30% in peripheral collisions |[T2l[T3l . provide a 
reasonable description of the data. 

Further hints on the 3/iff behaviour as a function of pi can be inferred from the centrality 
evolution of the (pj), extracted from a fit to the d 2 Nj/^/dydpj distributions. As presented in 
Fig j2](left), (pt) decreases towards central collisions, pointing to a smaller suppression of low pj 
J/tff. On the contrary, PHENIX results IfTTI [T5l [TBI show a different trend, with a (p T ) increasing 
with centrality, confirming the different J /iff behaviour, versus -\fs, already shown in FigJTJleft). 

Further insight on charmonium production in Pb-Pb collisions can be achieved by comparing 
the iff(2S) yield to that of the J/ift. Results are presented as a double ratio of the if/(2S) to J/i/f 
yields in Pb-Pb and in pp collisions as a function of centrality and in two pj classes (0 < pi < 3 
GeV/c and 3 < pj < 8 GeV/c). Signals are extracted with the aforementioned fitting procedure, 
keeping the ^(2S) mass position and width fixed relative to the i/iff ones, and then corrected for 
the corresponding Ax e. For 0-20% centrality, the very low signal over background ratio (S/B) 
prevents the extraction of the i^(2S) yield at low pj, while in the higher pj bin (S/B~0.01) only 
an upper limit can be evaluated. Since many of the systematic uncertainties cancel out in the 
double ratio, the main contribution to the systematic error is due the signal extraction, being of 
the order of ~ 15-60% depending on the kinematic bin. The pp reference has been evaluated at 
V^=7 TeV. Therefore, we have included a ~15% contribution to the systematic uncertainty to 
take into account a possible Vs-dependence G f the tf/(2S)/]/if/ ratio evaluated by comparing CDF 
IfTTI . LHCb HI and CMS U5] results. Double ratio results are shown in Fig|2] (right). The large 
statistics and systematic uncertainties preclude the drawing of strong conclusions on the i/'(2S) 
behavior. A significant enhancement in the double ratio for more central collisions is not visible 
in the ALICE data. 

We have presented ALICE results on the J/if/ Raa as a function of centrality and pj in Pb-Pb 
collisions at ys/wv=2.76 TeV, at forward rapidity. The Raa shows a clear reduction of the J/i/f 
yield, with a negligible centrality dependence and a strong pj dependence, especially in central 
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Figure 2: Left: the J/</> (pr) is shown as a function of Np art and is compared with results from the PHENIX experiment. 
ALICE statistical errors are shown as lines, while uncorrelated and partially correlated systematic uncertainties are 
shown, respectively, as boxes and brackets around the points. Right: the double ratio i^(2S)/J/i^ in Pb-Pb and pp is shown 
versus centrality and compared to CMS values |20|. Statistical errors are shown as lines, while systematic uncertainties 
are shown as filled boxes. The statistical and systematic uncertainties on the pp reference is shown as dotted lines. 



collisions. The ]/if/ behaviour presents different features with respect to the one previously ob- 
served by lower energy experiments. These features can be qualitatively described by theoretical 
models which include (re)combination as an additional mechanism for J/tf/ production. Further 
insight in the understanding ot the i/if/ behaviour still needs a precise knowledge of the cold 
nuclear matter effects, which will be studied in the incoming p-A data taking. 
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